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Abstract 

Superstring theory and other supersymmetric theories predict the existence of relatively hght, 
weakly interacting scalar particles, called moduli, with a universal form of coupling to matter. Such 
particles can be emitted from cusps of cosmic strings, where extremely large Lorentz factors are 
achieved momentarily. Highly boosted modulus bursts emanating from cusps subsequently decay 
into gluons, they generate parton cascades which in turn produce large numbers of pions and then 
neutrinos. Due to very large Lorentz factors, extremely high energy neutrinos, up to the Planck 
scale and above, are produced. For some model parameters, the predicted flux of neutrinos with 
energies > 10^^ eV is observable by JEM-EUSO and by the future large radio detectors LOFAR 
and SKA. 
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I. INTRODUCTION 



Cosmic strings could be formed as topological defects in the early universe. They are 
predicted in a wide class of particle physics models and can produce a variety of observational 
effects. These include gravitational lensing, linear discontinuities in the cosmic microwave 
background, and gravitational radiation, both in the form of a stochastic background and 
localized bursts. (For a review of cosmic strings, see, e.g., [H, H]-) 

Strings predicted in many grand unified models respond to external electromagnetic fields 
as thin superconducting wires As they move through cosmic magnetic fields, such strings 
develop electric currents. Oscillating loops of current-carrying string emit highly boosted 
charged particles from cusps - short segments where the string velocity momentarily gets 
very close to the speed of light. The emitted particles and their decay products can then be 
observed as high-energy cosmic rays j3] and gamma ray bursts jl-01- 

A phenomenon closely related to string superconductivity is the development of a bosonic 
condensate around the string core [sj . For example, a condensate of Standard Model Higgs 
particles could form around strings in some models. The Higgses would then be copiously 
produced at cusps, and their decay products could reach the Earth as cosmic rays |8|. 

Here we shall discuss an alternative mechanism of cosmic ray production, which does 
not assume string superconductivity or Higgs condensates. It relies on the existence of 
moduli - relatively light, weakly coupled scalar fields, predicted in supersymmetric particle 
theories, including string theory. Moduli would be copiously radiated by oscillating loops 
of string at early cosmic times, when the loops are smaller than the modulus Compton 
wavelenth, L < 1/m, and their frequency of oscillation is greater than the modulus mass. 
The emitted moduli may affect the big bang nucleosynthesis as they decay into photons 
and baryons, contribute to dark matter and to diffuse gamma ray background, resulting 
in stringent constraints on both the cosmic string tension and the modulus mass, when 
moduli are assumed to have gravitational-strength couplings to matter j9-12|. However, the 
couplings may in fact be much stronger, in which case the constraints from moduli radiation 
may be significantly relaxe d |l2| . Such strongly coupled moduli appear to be quite generic 
in string theory landscape |13l-ll8|. and this case is of particular interest for production of 
Extremely High Energy (EHE) cosmic rays and neutrinos. 

At later times, moduli can only be emitted from cusps, resulting in sharp bursts of high- 
energy moduli. Eventually moduli decay into standard model particles, and their decay 
products can be observed as cosmic rays with energies above lO^^ eV. 

The great interest of cosmic strings, and more generically topological defects, to high 
energy neutrino astronomy is based on tremendous energies of neutrinos accessible for these 
sources. While astrophysical sources can accelerate particles to energies 10^^ — 10^^ eV 
at most, topological defects can produce particles, including neutrinos, up to the Planck 
scale and above. Many observational methods of neutrino detection, in particular radio 
observations and observation of fluorescent light from space, are possible only above ~ 
10^° eV. Detection of such high energy neutrinos can by itself be considered as a signature 
of neutrinos from topological defects or other top-down scenarios. The production of EHE 
particles is a very generic property of topological defects, cosmic strings in particular, but 
large fluxes of such particles are produced only in exceptional cases 0], [lH - (22| . 

In this paper, we shall treat the modulus mass and coupling constant and the string ten- 
sion as free parameters. We shall estimate the EHE neutrino flux resulting from modulus 
decays and indicate some values of the parameters that can yield observable fluxes. The 
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paper is organized as follows. In Section II, we review modulus emission from cosmic string 
cusps (a more detailed derivation is given in the Appendix). In Section III, we discuss mod- 
ulus decay, EHE neutrino production, including beaming, and propagation in the universe. 
In Section IV, we review the size distribution of cosmic string loops and calculate the rate of 
bursts and the diffuse flux of EHE neutrinos. We also discuss here the upper bound on the 
neutrino flux, resulting from the diffuse gamma ray background observations. At the end of 
that section we give two illustrative examples of neutrino fluxes for different values of the 
model parameters. Finally, we discuss EHE neutrino detection. Conclusions are presented 
in Section V. 



II. MODULUS RADIATION FROM STRINGS 



The effective action for a modulus field interacting with a cosmic string of tension fi is 
given by (oj 



S 



d X 



i(V0)2 + lmV + 



47ra 



nir. 



0t; 



fi d cta/— 7, 



(1) 



where 7 is the determinant of the induced worldsheet metric ■jab = 9tJLuX%X% -^^(c"; ''") is 
the string worldsheet, is the trace of the energy momentum tensor of the string, a is the 
modulus coupling constant, m is the modulus mass and mp is the Planck mass. For a ~ 1, 
the modulus coupling to matter is suppressed by the Planck scale. Here, we treat a as a 
free parameter and are mainly interested in a ^ 1. Then, the mass scale characterizing the 
modulus interactions is ~ nip/ a <^ nip. Values as large as a ~ 10^^ have been discussed in 
the literature 17 . 



The modulus field equation has the form 

(V'-m2)0(x) 

with 



At: a 



nir, 



(2) 



T;(x) = -2/i / drdav^'^^(a;" -a;"(a,r)). 



(3) 

The power spectrum of modulus radiation from an oscillating loop of string can be de- 
composed in Fourier modes as [9| 



dPr, Ga" 



■u;„A;|T(k,u;„)| , 



dn 2tt 

where G is the Newton's constant, ujn = ^/W^^-w? = Ann/L, L is the length of the loop 



T(k, a;„) = — — d x dadr^ —^5 (x" — x"(cr, r))e 



(4) 



(5) 



and = (a;„, k). 

We shall be interested in the modulus emission from large loops of string, having length 
L ^ m~^. In this case, the characteristic frequency of loop oscillation is a; ~ 1/L -C m, 
so modulus production is suppressed, except in the vicinity of cusps, where extremely high 



3 



frequencies can be reached in a locahzed portion of the loop for a brief period of time. 
Lorentz factors greater than^ 7 are reached in a fraction of the loop of invariant length 
AL ~ L/7. 

The spectrum of resulting particle bursts can be found by expanding X^((T, r) near a 
cusp [sl, [lo|.^ One finds that the number of moduli emitted in a single burst with momenta 
k in the interval dk (in the center of mass frame of the loop) is given by 

dN{k) ~ a^Gfi^L^/^k-''/^dk. (6) 

This distribution applies for k > k^, where 



kc ~ —mVniL. (7) 

At smaller k the distribution is strongly suppressed, dN ^ 0. 

The dominant contribution to the modulus emission comes from the lower momentum 
cutoff kmin ~ kc, SO the total number of moduli per burst is 

N ~ (8) 

The particles come from a portion of the loop that reaches Lorentz factors in excess of 

7c ~ kc/m ~ ^VmL, (9) 

and are emitted into a narrow opening angle -t^c around the direction of the string velocity 
V at the cusp, 

t9e ~ Ic^ ~ 4(mL)-^/l (10) 
The total power of modulus radiation can be similarly calculated as 

Pm ~ a^Gfi'L-'/'k:'/^ ~ (11) 

y/mL 

The loops also radiate gravitational waves with the power 

Pg ~ rGfl', (12) 

where F ^ 50 [l|. Pg ^ Pm when L L^, which is given by 

~ T-^a^m-\ (13) 
The lifetime of a loop which mainly radiates gravitationally is 



^ From here on, and until Appendix, we use the notation 7 only for the Lorentz factor. 

^ A detailed analysis of particle emission from cusps is given in the Appendix, confirming the results obtained 

in 13 0. 
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which imphes that the characteristic size of the smallest (and most numerous) loops surviving 
at time t is 

L^, ~ TGfit. (15) 

On the other hand, modulus radiation dominates when Pg < and the loop lifetime is 
given by 



Pm a'^Gjj 
The corresponding minimum loop size at time t is 

LZn - a^/^(G/.)^/^m-V3^2/3_ (17) 
The transition between the two regimes occurs at 

Therefore, the minimum loop length is given by f lT5]) for t > and by f[T71) for t < t^. The 
redshift corresponding to is given by 

~ ^2a-«/3(G^)2/3(mto)'/^ (19) 

numerically 

z, ~ 4007715/^07 (20) 

where 

7775 = 777/10^^6^, aj = a/10\ /i_2o = G/i/lO"^", (21) 
and the fiducial values have been chosen anticipating the results in Sec. IV. 



III. NEUTRINO PRODUCTION AND PROPAGATION 



In this section we address some problems in neutrino physics relevant for future con- 
sideration, namely, the neutrino horizon, the energy spectrum of neutrinos produced by a 
modulus decay, the boost of this spectrum by the cusp Lorentz factor, and some others. 

We start with a note about accuracy of our calculations. 

The main purpose of our work is a discussion of the principle features of the phenomenon, 
i.e., EHE neutrino production by moduli from cosmic strings, not an accurate numerical 
evaluation of the neutrino fluxes and their detection rates. In particular, our aim is to express 
the results in the form of analytical formulae, so that the dependence on input parameters 
can be easily seen. For this purpose we make the following simplifying assumptions. 

We use the CDM cosmo logical model with A = and flm + fir = ^ and use Hq = 
72 km/sMpc, to = 4.3 x 10^^ s, teq = 2.4 x lO^^ i ^ ^ 3200, the scale factor in the 
radiation and matter dominated eras are ar{t) oc t^^^ and am{t) oc t"^^^. The corresponding 
time-redshift relations are respectively given by (t/to) = (1 + ZeqY^'^{l + z)^"^ and (t/to) = 

(1 + z)-3/2. 

For convenience of calculations we assume that at the decay of a modulus at rest the 
neutrino spectrum is oc while in reality this spectrum is not power-law and is approx- 
imately proportional to E~^''^ only for a very large mass of the decaying particle. 
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A. Neutrino horizon 



As they propagate through the universe, the UHE neutrinos Vi (or antineutrino Vi) with 
i = e, iJ,,T are absorbed or loose energy in the following three reactions: 

(i) Vi + fi ^ Qa + Qa, {a) Ui + Ui I + I, {in) z/j + z/j -)■ i^j + Uj, (22) 

where q = Ua, da, Sa, Ca, ha are quarks with a = 1, 2, 3 being color indices and I = e, fi, t 
are lepton flavors. Reactions {i) include only s-channel, and {ii) may include both s and t 
channels. For a rough estimate we can use the following approximation for the cross-section, 

{(N/Tr)G'is at s < rnL 
..r/ ^r^2 2 r ' (23) 

{N/7r)Gj,m^ at s > 

where Gp = 1-17 x 10~^ GeV~^ is the Fermi constant, s{z) = 2Eymy{l + z) is the center 
of mass energy squared at redshift z, Ey is the neutrino energy at the present epoch, m,^ ~ 
0.1 — 0.2 eV is the assumed neutrino mass, and iV ~ 10 — 15. 

UHE neutrinos are absorbed or loose energy in collisions with relic background neutrinos 
whose space number density is = 56(l+z)^ cm~^ and kinetic energy is e,, = 3.15T{l+z) = 
5.29 X 10-^(1 + z) eV. 

Neutrino horizon , i.e., the maximum redshift z^, from which neutrino with the observed 
energy Ei, can arrive, is calculated as 

/ dt a,{z) n,{z) = 1, (24) 

Jz=0 

where 

dt = io(l + zr^^^dz. (25) 

At the highest neutrino energies, when a ~ amax ~ {N/n)Gpmly , z^ ~ 1.5 x 10^. At 
energies below 2 x 10^^ GeV, 

Zy ~ 2.5 X 10^(^/10^^ GeV)-2/^ (26) 

For energies of interest in this paper we shall use z^ ~ 200 at E^ > 10^° eV. Detectable 
UHE and EHE neutrinos are produced in the matter-dominated epoch. 

It is interesting to note that the modulus-string model allows to probe the earliest universe 
with the help of EHE neutrinos: e.g., for superconducting strings j3] the maximum redshift 

is ^max ~ 3. 

B. Modulus decay and neutrino spectrum 

The rate of decay of a modulus into the Standard Model (SM) gauge bosons can be 
estimated as ^ 

Fo-ricjA/f— ) m^ (27) 
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where hsm = 12 is the total number of spin degrees of freedom for all SM gauge bosons, m 



is the modulus mass, and we assume interaction of the form [14 



A„t <PF,,F>^\ (28) 

The mean lifetime of the modulus in its rest frame is then 

To ~ 8.1 X 10"^^m^^af 2 s. (29) 

For a wide range of parameters and a, the lifetimes of moduli are short even after a 
large Lorentz boost. In our main scenario, the neutrino-producing moduli are born within 
the neutrino horizon and decay almost momentarily there. 

However, in principle the redshifts z > Zi, can also contribute to the neutrino flux at 
z = 0, if the lifetime of boosted moduli is long and they can decay a.t z < z^. In the analysis 
below we argue that such a scenario is disfavored. 

The lifetimes of moduli emitted from cusps are boosted by large Lorentz factors. A 
modulus emitted with a Lorentz factor 70 = fc/m at redshift z and decaying at redshift Zd 
has a lifetime ^ 

t{z) ~ ro7o— — ^. (30) 

In order for neutrinos to reach the Earth, they should be produced within the neutrino 
horizon at redshifts z^ < z^. Moduli emitted from cusps at z > Zj, can therefore yield 
observable events only if they have large enough lifetime, allowing them to survive until 
they reach z^,. This gives the condition 

riz) ~ ro7o^^ > t{z^) ^ to(l + z^)-'^', (31) 

where in the last step we used the fact that in the energy range of interest Zi, < Zeq. 

For z > z,y and using > 1 and a? > 1 (which is necessary for detectable neutrino flux; 
see Sec. IV), we obtain 

> to/To 3.8 X 10^^ 

which is too large at any Zd- 

In what follows we shall consider only cusp events occurring at z < z^ < z^q- There are 
no restrictions on the modulus lifetime in this case, except that it should be short enough for 
sufficient fraction of moduli to decay before they reach the Earth. This is always satisfied 
in the parameter range of interest. 

The decay channel relevant for the neutrino production is the decay into gluons, via the 
modulus-gluon interaction of the form fl28l) . The primary gluons from the modulus decay 
initiate the quark-gluon cascade, which turns into hadrons, mostly in pions, and then to 
neutrinos. 

To simplify calculations and analysis, we shall assume the neutrino production spectrum 
dN/dE oc -E"^, close to the power-law approximation E~^'^ obtained for a large mass of the 
decaying particle using Monte Carlo simulation and DGLAP method (23 |. 

The neutrino spectrum from a modulus at rest is then 

= i^{E,,m) ^ ^fr^K^, (33) 
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where m is modulus mass, is neutrino energy, b^, is given in terms of the ratio of maximum 
^max ^ 0.1m and minimum e™™ neutrino energy, 

K = [HeTVCn]'\ (34) 

and /tt ~ 1 and 1/2 are the fractions of energy transferred from the modulus to pions and 
from pions to neutrinos, respectively. 

The spectrum of low-energy neutrinos is a model-dependent feature, but generically sup- 
pression of this spectrum is provided by suppression of soft gluon emission due to coherence 
effect in the parton cascade. It results in the Gaussian peak in the spectrum of pions, parents 
of neutrinos. We describe effects of low-energy suppression of neutrino spectrum introducing 
formally the energy e™™ in Eq.( 1M|) . where the suppression starts, and refer to it as "the 
minimal energy" . Determining the value of e™" would require numerical calculations. Here, 
we shall parametrize 

£ = 71 GeV. (35) 



C. Lorentz boost and beaming 

Emerging from a cusp segment, a modulus obtains very large Lorentz factor 7 corre- 
sponding to the point of exit. The typical Lorentz factor, as it is calculated below, reaches 
7 ~ 10^^ — 10^'^. The neutrino energy in the modulus rest system is boosted as 

E = 7^,(1 + /3 cos ^,), (36) 

where /3 = v/c, and "(9* is the angle between the directions of neutrino motion and of the 
boost in the rest frame of the modulus. 

First of all we calculate how the neutrino spectrum f l33|) changes under the transformation 
f p6|l . For this we use the conservation of the number of particles before and after the Lorentz 
boost: 

dN^{E,,^,) = h^m^^ = dN{E,E^). (37) 

Using 

dE 

dfl^, = 2n (icosi?* = 2n , (38) 

IpE^ 



we obtain in terms of the new variables E and E^, : 

d'^N , m 



dEdE, *A-i(5El' 
Integrating with respect to dE^ we have 



(39) 



dN _ m 1 



where the minimum neutrino energy for a fixed E is 



E 

-'-'* mm 



E/[j{l + f3)] if i?>7(l + /3)C- 

(^) = <^ . . , (41) 

[eT^ if ^<7(1 + /3)C 
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Finally, we obtain for /3 ^ 1 and for decay at epoch z 



if the neutrino energy at the present epoch is E > 7(1 + For E < 7(1 + 

dN/dE does not depend on the neutrino energy E : 

dN ^ . ^ . 1 k 1 , , 

— ^Ui?,A;,.) = -^-^. (43) 

In both formulae above, k = 7m is the modulus energy. 

We now briefly discuss the effect of beaming. Strong beaming of the produced particles 
is a remarkable feature of the cusp models, which provides interesting observational conse- 
quences. For moduli with the range of masses considered here the beaming is very strong 
(see Section [IT])- When gravitational radiation dominates, the Lorentz factor aX z = can 
be estimated as 7c ~ 4.5 x 10^^ 'f^H'^ 1^^-20 ■ When modulus radiation dominates, the Lorentz 
factor at z^ is 7c = Q;^/4r = 5 x lO^^cty. As a typical Lorentz factor for neutrino production 

one may consider that at neutrino horizon Zu ~ 200, 7c = 8.5 x lO^^mg'^^/i^^o- 

In the frame where the modulus is at rest, neutrinos are emitted isotropically. After a 
Lorentz boost, most of them move within a narrow cone with ~ I/7 in the direction of 
the boost, with energies E ~ 7-E*. However, neutrinos which are emitted in the rest frame 
within a narrow cone with -(9'^^ = I/7 in the backward direction are moving after the boost 
within a wide angle ~ 7r/2 in the direction opposite to the narrow high-energy jet and 
typically have very low energies E ~ E^/'-f. The total number of these neutrinos is 47^ times 
smaller than that in the high-energy jet, and they are undetectable because of their small 
number and low energies. The typical neutrino energies in high-energy beam is -E ~ lE^,, 
but low-energy neutrinos with E ~ E^/'y are also present there. 

Formally, the minimum neutrino energy is extremely low, E ~ £:™"/7, but the number 
of such neutrinos is very small. However, when neutrino spectrum is changing from at 
high energies to dN/dE oc const, as in Eq. ( H3|l . the detectability of neutrinos is sharply 
decreased, and thus -Emin ~ 7(^)er™ can be regarded as an effective low-energy end of the 
spectrum at epoch z. To estimate the low-energy spectrum cutoff for neutrinos generated 
at epoch z^, ~ 200 and observed now, we use the parametrization fl35p and our estimate for 
the Lorentz factor at z^, ~ 200, 7c = 8.5 x lO^^rng'^^yui^V Taking into account the redshift of 
the neutrino energy, we find -Emin ~ 4.3 x lO^em^^^/i^^^Q GeV. This estimate demonstrates 
the remarkable feature of our model: the predominant generation of the highest energy 
neutrinos. A more realistic estimate of -Emin for the diffuse neutrino flux will be made in 
Section IVCl 

EHE neutrinos propagate as a jet in a cone with an opening angle ^9 ~ I/7. The duration 
of the neutrino pulse is very short, r ~ L/'y^, and all neutrinos reach the detector almost 
simultaneously, due to the smallness of the neutrino mass. The effective area illuminated by 
arriving neutrinos is much larger than the area controlled by the detector, but simultaneous 
appearance of two-three showers in the field of view of a large detector is possible for some 
parameter values (see 0] for discussion and calculations). 

Beaming is a property of all particles emitted from a cusp, in particular gamma-rays. In 
some models part of observed GRBs may be produced by cosmic strings ^] . Since the total 
number of strings in the Milky Way is tremendously large {N ~ 10^), one may expect that 
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UHE GRBs from the Milky Way may be observable. In fact, the rate of predicted bursts 
is too strongly suppressed by the beaming factor f2/47r = 1/(47^) to be detectable. The 
backward component is distributed within a wide sohd angle and is not suppressed by the 
this factor, but the total number of photons and their energies are too small for detection. 



IV. NEUTRINO BURSTS FROM MODULI 
A. Loop Distribution 

The predicted flux of EHE neutrinos depends on the typical length of loops produced by 
the string network. The typical loop sizes have been a subject of much recent debate, with 
different simulations 25-3^ and analytic studies [s^, S^] yielding different answers. Here we 
shall adopt the picture suggested by the largest, and in our view, the most accurate string 
simulations performed to date |32l]. According to this picture, the characteristic length of 
loops formed at cosmic time t is given by the scaling relation 

L ~ f3t, (44) 

with /3 ~ 0.1. 

The number density of loops with lengths in the interval from L to L+dL can be expressed 
as n{L,t)dL. Of greatest interest to us are the loops that formed during the radiation era 
{t < teq) and still survive at t > teq. The density of such loops is given by 

n{L, t)dL ~ p-\{(3t,gy/H-^L-^/^dL, (45) 

where p is the string reconnection probability and C ~ 16 is the parameter characterizing 
the density of infinite strings with p=l, poo = C/^/^^- 

The dependence of the loop density on p is somewhat uncertain and can only be de- 
termined by large-scale numerical simulations. Here we have adopted the p~^ dependence 



suggested by simple arguments in, e.g., [35|, [361]. The reconnection probability is p = 1 



for ordinary cosmic strings. Its value for F- and D-strings of superstring theory has been 
estimated as ji^ 

10-3 <p<l. (46) 

The distribution f l45p applies for L in the range from the minimum length L^m to L^ax ~ 
(3teg. The lower cutoff L^m depends on whether the energy dissipation of loops is dominated 
by gravitational or by modulus radiation. It is given by f[T^ for z < z^, and by f[T7l) for 
z^ < z < Zeq, with z^ from Eq. (l20!) . For > z^q, the dominant energy loss is gravitational 
radiation and Eq. (fT5|) for Lmin applies in the entire range of interest. 

The string motion is overdamped at early cosmic times, as a result of friction due to 
particle scattering on moving strings. The overdamped epoch ends at [l| 

tdamp ~ {G^lY'^tp, (47) 

where tp is the Planck time. In the above analysis we have assumed that loops of interest 
to us are formed at t > tdamp- The corresponding condition is 

Lmin{t) ^ Ptdamp- (48) 
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This bound assumes that the strings have non-neghgible interactions with the standard 
model particles, so it may not apply to F- or D-strings of superstring theory. In any case, 
we have verified that (148|) is satisfied for parameter values that give a detectable fiux of 
neutrinos. 

B. Gravity- and moduli-dominated regimes. Restrictions imposed by 

The value of 2r* given by Eq. fl20l) marks the boundary between two regimes of string 
evolution. In the first regime, at z < z^, the string energy loss is dominated by gravitational 
radiation. This includes the entire relevant range of z for z^, > Zy and the range < z < 2^, 
for z^, < Zy. We shall call it the gravity- dominated regime. 

The second regime is dominated by the modulus radiation; we shall call it the moduli- 
dominated regime. It corresponds to the redshift interval z^ < z < and exists only when 

•2^* Zfy . 

It is often convenient to perform the calculations for the fixed value of z^, considering 
it as a free parameter. In this case the space of three physical parameters a7,m5 and fi is 
restricted by Eq. (120|) as 400^7 ^'^^rng'^^/x^^^Q = z^,, which we will use in the form 

/i-20 = (^*/400)3/2a4 ^-1. (49) 

Then our calculated quantities, such as neutrino fiuxes and characteristic energies, will 
depend on two parameters, a? and m^, and the fixed value of z^. 

The value of fixed z^, characterizes "the model" . One should distinguish two major classes 
of models: with high z^ > z^ , which corresponds to gravity-dominated regime, and with 
low z^ < Zy, which includes both the gravity-dominated regime at < z < z*, and the 
moduli-dominated regime at z^ < z < Zy. 

C. Neutrino Flux 

As we argued in Section IIII Bt the neutrino-producing moduli are born and decay within 
the neutrino horizon at z < z^. We shall first consider the high models with 2;* > Zy 
(gravity-dominated regime) and then study the case z^, < z^, which includes both the moduli- 
dominated regime at z^ < z < z^ and the gravity-dominated regime at < z < z^,. 

The neutrino fiux can be most generally calculated as 

^"(^- ^' = i / TTT '"'^^"^ I (uk) ^-f^- 

where the proper volume for the matter-dominated epoch is 

dV{z) = 5A7rtl[{l + zfl'' - \f{\ + zY^^I'^dz , (51) 

the rate of bursts is dNh = n{L)dL/{L/2) with n{L)dL being the number density of loops 
with length L in the interval dL (see Eq. (HSj)), the number of moduli dNx emitted in a 
burst with momenta k in the interval dk is given by Eq.([H]), ^^A:/47r is the probability that 
a randomly oriented burst is directed to the observer, fifcr^(z) is the area of the irradiated 
spot at the observer's location, and ^^{E, z, k) is the spectrum of neutrinos from decay of a 
modulus with momentum /c, given by Eqs. (H2|) and (H3|) . 
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Integration in Eq. (l50i) goes over k, L and z. For integration over k and L only lower 
limits are essential, and they are given by fcmin = my/mL and Lmin from Eq. f lTSj) or Eq. f|T7|) 
for the gravity-dominated and moduli-dominated regimes, respectively. In the case z^ > Zi,, 
we have ^max = z^, while 2;min is determined by the rate of bursts or by minimum energy of 
neutrino -Emin(^) at epoch z as explained below. 

Consider first the limit z^i^ imposed by the rate of bursts. The average rate of bursts 
Nfy{< z) that occur in the redshift interval between and 2; is a growing function of z, and 
we can define Zb as the redshift at which this rate is a few bursts per year. No bursts will be 
detected from z ^ Zb, so we should introduce a lower cutoff of z-integration at Zmin = Zb- 
The rate of burst Nb is calculated in Section HV DI and Zb is found to be small in the parameter 
range that we are considering here. Hence, the condition z > Zb does not yield a significant 
constraint for the 2;- integration in Eq. (l50ll . 

Another constraint to consider is that the energy E of neutrinos should be above the 
minimal energy, E > E^in{z) = 7c(-z)£^™"/(l + z), where 7c(-2) = |["^-Lmin(-2)]^''^ is the 
characteristic Lorentz factor for loops of minimal length L^ia{z) at epoch z. For the gravity- 
dominated regime we are considering here 

7c(^) = 4.5 X lO^^^mf /.^/^(l + (52) 

To proceed, it will be convenient to use z^ in place of the string tension as a free parameter. 
With the aid of Eq. (H9|) we have 

E^U^) = E,a^,{zjz,f/\l + zr"\ (53) 

where 

£"0 = 2.7 X lO^^e GeV, (54) 

and e is the parametrization factor introduced in Eq. fl35l) . The lower bound of 2;- integration 
Zmin can now be found as the value of z for which Emin{z) = E, 

1 + ZrmniE) = {E^j E^' {z^j z^f. (55) 
Integrating Eq. (!50|) over z from Zmin{E) to z^, since we assumed above z^, > z^, we obtain 

The calculated flux for "normalizing" set of parameters p = 1 (ordinary strings) and 
ar = ms = z^:/zi, = 1 is shown in Fig. [Uby curve "theor.3". This flux is low and detectable 
only by SKA. The largest flux in Fig. [T]is presented by curve "theor.l" for the parameters 
p = 1, ai = 2, ms = 0.1, z^,/z^ = 1 . It is close to upper limit shown by curve 
cascade", and is detectable by JEM-EUSO, LOFAR and SKA. Here and everywhere below 
we assume e = 1. 

The maximum energy of neutrinos is determined by -Emax ~ 27e™'^^ at generation and 
can be extremely large, but the flux of these neutrinos is suppressed as E""^. 

We shall now consider the low 2* models with z^, < z^, and calculate flrst the neutrino flux 
generated in the redshift interval z.^ < z < Zy, where energy losses are moduli-dominated. 
Then we calculate flux in the interval < 2; < z*, where the gravitational radiation domi- 
nates. 



1 ~l~ Zffiin (£') 
1 + 



1/2 



GeV/cm^s sr. (56) 
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For the interval < z < and fixed z*, the neutrino flux is calculated using Eq. (|50 
and the parameter restriction in the form yU%^o ~ (2;*/400)ay^^m5 As a result we have 



with = 1.8 X 10^1° GeV/cm^s sr. 

The lower limit of integration in Eq. (J571) is obtained as above from Ejam{z) 
7c(-2)£^min/(l + z). Usiug the condition E > E^i^{z) and Zmin > ^* one obtains 

(Z^ at E > Emin 

(E/Eo)"'/'"?^'^^ at E^^<E<E^^, 
z„ sX E < E^^„ 



(57) 



(58) 



with Ernin = 1-3 X 10^^ea^{z^/50y/'^ eV, = 1 x lO^^ea^jidO/z.y eV, and Eq = 5 x 

lO^^e eV. Finally we have 



1 + ^^ 

1 ~)~ ^min 



1/2 



- 1 



(59) 



Using in Eq. ( 15^ 2;jnin(-£') = z^ at E > Ernin (see Eq. (1^ ) one finds the high energy (HE) 
asymptotic of the fiux in the moduli-dominated regime 



E\J^{E) = 1.8 X 10-^%-^a^^m^^^\l + z 



a/2 



1 + 



1/2 



GeV cm"^ s"^ sr-\ (60) 



valid at _E > -Emin- 



We have to add also the neutrino fiux generated in the interval < z < where 
the gravitational energy losses dominate. In this case Lmin is given by Eq. ([T^ . For HE 
asymptotic one finds 



E^J^{E) = 1.8 X 10->-^a2^^^/^(l + z^y/^ GeV cm'^ s 



(61) 



For an easily understandable reason, the HE asymptotic here is {z^^/z^y^'^ lower than in 
the case > Zy given by Eq. f lSB]) . Less trivial is the coincidence of HE asymptotic in the 
gravitational regime (0, z^,) given by Eq. f ET]) and in the moduli-dominated regime (z*, Zi,) 
given by Eq. ( 160|) . It is explained by the fact that dominant contributions in both cases 
are given by epochs with redshift 2;*. The general conclusion about neutrino production is 
therefore increasing the fiux with growth of z^, until it reaches z^. Thus, the high z^ models 
predict the largest neutrino fiuxes. 

Apart from this, more efficient neutrino detection in high models, follows from a lower 
cutoff energy -Emin- Indeed, the minimal energy of neutrinos for both regimes (0. z*) and 
{z^, Zy) is determined by the same expression 



EjaixiiZif'] 



lc{z*) 
1 + Z^ 



la^ el 



4 11 + ^/ 



(62) 
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while for the gravity-dominated regime with > Zy it is determined by Zy as 

E^^{zy) = ^^er- (63) 

Increasing z^, in Eq. ( l62l) we decrease -Emm, which is favorable for detection by JEM-EUSO. 

One may also see that moduli-dominated regime gives subdominant effect as compared 
with gravity-dominated 



D. The Rate of Bursts 

The rate of bursts is not a physically measured quantity, but it can serve as indicator of 
detectability of the burst radiation. 

The rate of cusp bursts with their cones of radiation directed to the observer is given by 

^ /-^n^i^ £ (64) 

J 1 + z L/2 Att' ^ ^ 

where dV{z) and n{L)dL are the same as in Eq. (1501) . The quantity 0/47r gives the proba- 
bility for the observer to be located within the cone of the cusp radiation. Such a location 
does not guarantee detection of this radiation, because the area of irradiated spot Qr"^ is 
much bigger than the size of the Solar system, but too low rate of bursts means that the 
signal is undetectable. To calculate ~ Tcd'^ we use ~ {kcL)~^/^ (see the discussion above 
Eq. (11031) in the Appendix for details). 

We calculate the rate of bursts for z < z^, when gravitational radiation dominates. The 
rate integrated from 2; = up to redshift z is given by 

M< z) = ^2^/3 p-\py- 10^^ -L /(,) = 1.1 X 10^ f^Zli' /(.) yr-\ (65) 
where 

i{z) = r dz'[{i + z'y/^ - i]2(i + z'y/\ (66) 

Jo 

Numerical values of I(z) are shown in Table I. 



TABLE I: Numerical values of integral I{z) 



z 


0.1 


0.5 


1.0 


2.0 


I(z) 


9.12 X 10"^ 


0.015 


0.165 


1.99 



In Section IIVCI Zmm = is defined from condition iV';,(< z^j) is a few per year. From 
Table [T] and Eq. (165|1 one can see that the rate of bursts is large enough even at redshift as 
small as 0.1, and thus the condition z > z^ does not impose a significant constraint for the 
integration of Eq. (1501) over z. 
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E. Cascade Upper Bound 



An upper bound on the neutrino flux follows from the observed diffuse HE gamma-ray 
background, since neutrino production via pion/kaon decays is accompanied by high energy 
electron and photon production. These electrons and photons, interacting with CMB and 
Extragalactic Background Light (EEL) photons, produce an electromagnetic cascade, whose 
energy density Ucas must not exceed that of the observed diffuse gamma radiation. This 



results in the upper limit on the diffuse neutrino flux j47|. With the assumption of E 



generation spectrum of neutrinos the cascade upper limit can be written as 



„ , .max 



An ln(Emax/^min) 

where i^max and £'min are the maximum and minimum neutrino energies, respectively, and 
oj^^ is the maximum cascade energy density allowed by observation of the isotropic diffuse 
gamma-radiation. 

According to recent Fermi-LAT observations |48l| this energy density is = 5.8 x 

10~^eV/cm^, as it follows from the analysis jiol, |50|. This limit results from a comparison 
of the cascade spectrum at 100 GeV with the Fermi data, while the lower energies give a 
weaker limit. For our case we assume that cosmological epochs with redshifts z > Zcas, 
when gamma-rays with > 30 GeV are absorbed, make a negligible contribution to the 
obtained upper limit. One can estimate Zcas in the following way. At z = a very sharp 
absorption of gamma-rays on CMB occurs at i?"''^(0) ~ 100 TeV. At epoch z this energy 
is {1 + z) lower. Taking into account the redshift of these photons, we estimate z^as as 
z,as ~ {E^^'/E.yl^ ~ 60. 

The energy density for the electromagnetic cascade radiation resulting from modulus 
decays can be expressed as 

where /^r ~ 1 and 1/2 are the fractions of energy transferred from the modulus to pions and 
from pions to electrons and photons, respectively, dt is given by Eq. (!25|l . the density of the 
loops n{L) is given by Eq. (H5|) . and the number of moduli emitted per burst N{k) is given 
by Eq. ([6]). The limits of integration in Eq. (!68|) are as follows : z^^i^ = 0, 2;max = Zcas ~ 60, 
kmin = {l/4:)mViTiL, and Lmin is given by Eqs. ( !T5|) and f|T7|) for the cases > Zcas and 
< Zcas, respectively. 

Consider first the case > Zcas, when gravitational radiation dominates. Performing the 
integration we obtain 

Ucas = 9 X 2-l/V^C/3^/'r-^M^a2^(l + ZcasY^', 

= 5.8 X 10-V^a?"^^^^^(-2cas/60)i/2 eV/cml (69) 

For z^: < Zcas we have to integrate over the interval < z < z^, where gravitational 
radiation dominates, and over interval z^, < z < Zcas, where moduli radiation prevails. For 
the first case (gravity-dominated regime at < z < z^ ) one obtains 

ojcas = 4.7 X IQ-^p-^a'^^ m^^^\z,/AOy/^ eV/cml (70) 
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For the second case (moduli-dominated regime at < z < Zcas) we have using Eq. ( H9|) 

Ucas = 4.7 X m-'^\z, /A0y/^[1 - (z,/ z^as)'^^] eV/cml (71) 

Note that this is the same as (170 p . apart from the last factor. The reason is that in both cases 
the main contribution to z-integration comes from 2; ~ z*, like in case already discussed for 
neutrino fluxes. 

For reasonable values of z*, all three Ucas given by f l69|) - flTTj) are less than maximally 
allowed u^^f" = 5.8 x 10"'^ eV/cm^. 

Two remarks about the cascade limit for our model are now in order. 

This limit is not strictly enforced because cascade production occurs at very large red- 
shifts, while Ucas is constrained mainly due to the highest energy cascade photons with 

^ 100 GeV (see Fig. 1 in j49|). These photons are absorbed by EBL at earlier cosmo- 
logical epochs, and therefore in cosmic string models, where the main part of the cascade 
energy density is produced at large redshifts, the constraint on cOcas is weaker and values 
higher than 5.8 x 10"'' eV/cm^ are allowed. Therefore, neutrino fluxes above -E'~^-cascade 
upper limit taken from [49| and shown in Fig. [1] are not necessarily excluded. 

One may also expect that beaming of cascades may reduce the efficiency of the cascade 
limit. The calculation of beam widening in the ambient magnetic field shows that narrow 
beaming survives only for very weak magnetic fields of order 10~^^ G. 



F. Superstring-motivated example for high z^, model with > z^, 



Because of the simplifying assumptions adopted in this paper, we cannot reliably deter- 
mine the domain in the parameter space of a, m and G/i, which yields detectable fluxes of 
EHE neutrinos. 

Instead, we shall discuss some illustrative examples of parameter choices, for which neu- 
trino detection in future experiments is possible. Among such projects we shall consider 
detectors aimed at the highest energy neutrinos above 10^° eV: JEM-EUSO 37|, Anita II 
m, LOFAR lli, |39|, SKA|38|, |io|, and LORD |4l|. JEM-EUSO is based on space detec- 



tion of fluorescent light from the showers in the atmosphere. All others are based on radio 
detection of the showers due to the Askarian effect |42|. The sensitivities of three of these 
experiments, JEM-EUSO (expected to be launched in 2015), LOFAR and SKA, are shown 
in Fig. [TJ 

First we shall consider a superstring-motivated example for cosmic strings evolving in 
the gravity-dominated regime, assuming that z^ > Zy. We flx plausible parameters, both for 
strings and for particle physics. Speciflcally, we consider the large volume string compacti- 
flcation model 13||, which is characterized by an intermediate string scale iris ~ 10^^ GeV 
and a TeV-scale supersymmetry (SUSY) breaking. The hierarchy between the Planck and 
SUSY breaking scales in this model is due to an exponentially large volume of the compact 
extra dimensions, Vcomp = V/f , where V ~ 10^^ and Is ~ mj^ is the string length scale. 

Apart from the volume modulus, which has gravitational strength couplings to ordinary 
matter, the other Kahler moduli have large couplings of the order ^ 



a ~ 



V. 



(72) 



With V ~ 10^^ we have a ~ 10^'^ 
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The modulus masses are given by [l 



m ~ -^rup. (73) 

It is useful to parametrize m in terms of a. Since we are interested in a ~ 3 x 10^, the factor 
InV can be replaced by 35. Hence, we obtain for the mass 

m ~ 35mpa~^. (74) 

For a ~ 3 X 10^ we have m ~ 4 x 10^ GeV, and we fix /i_2o ~ 10, to ensure > z^. The last 
condition means choosing the basic cosmic string parameter, the symmetry breaking string 
scale ri = ^JJi ~ 3.9 x 10^ GeV, i.e. G/i ~ 1 x 10~^^. With this choice of parameters, z^, from 
Eq. (1201) is given by 2;* ~ 250, i.e it is larger than z^^^ ~ 60 and than the neutrino horizon 
Zy ~ 200. For this case, the neutrino flux is given by Eq. and the cascade energy 
density by Eq. fl69|) . With the parameters ay and as indicated above, the neutrino 
flux is E'^J{E) = 1.3 X 10~^p~^ GeV/cm^ s sr and the cascade energy density is Ucas ~ 
2.6 X lO'V^ eV/cm^. 

This neutrino flux is shown in Fig. [T]by the curve "theor.2". It is detectable by LOFAR 
and SKA, but not by JEM-EUSO. 

We can modify this model choosing parameters providing a larger neutrino flux. For this 
we flx ^7 = 2, ms = 0.1 and keep z^: = 250, i.e. /i_2o ~ 80 and r] ^ Ix 10^° GeV. The cascade 
energy density, Wcas = 7.3 x 10~^p~^ eV/cm^, is safely below u^"^. This model satisfies the 
restrictions obtained in il2!|, and in fact we can further increase the flux by decreasing 777.5. At 
high energy limit the calculated neutrino flux is given by E'^Jy{E) ^ 3.2 x 10~^p~^ GeV/cm^ 
s sr. This flux is shown in Fig. [Uby curve theor.l, it is detectable by JEM-EUSO, LOFAR 
and SKA. 



G. An illustrative example for low model with < Zjy 

Let us now assume that < z^y. Then, in the redshift interval z^ < z < z„ the modulus 
radiation dominates. In the remaining interval < z < z* the gravitational energy losses 
are dominant. The total neutrino flux is given by the sum of fluxes generated from both 
intervals. A value of z^, flxes a model. Here we try to flnd a model with a detectable neutrino 
flux. 

The neutrino flux from the moduli-dominated interval {z^, ,Zi,), in the HE asymptotic 
regime, is given by Eq. f lBU]) . This equation becomes valid just above the low-energy steep- 
ening, i.e. at 

E > E^in = 1 X 10^*^ a^j {50/z,f eV. (75) 

For the flux from the gravity-dominated interval (0, z^,), the HE asymptotic is given by 
Eq. ( 16T|) . with a low-energy cutoff approximately at the same energy as above. Summing up 
both components, we obtain 



E'Jy{E) = 1.8 X 10-> 



1/2, 



1 + 2;* 



\l/2 



1 + ^* 



1/2- 



GeV cm^" s"' sr^^ (76) 



To maximize the flux without strongly increasing the low-energy steepening threshold 
-E'min given by Eq. ( 175|) . one can use the parameters: = 2, 7ri5 = 0.1 and z^ = 100. As a 
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result we obtain .Emm = 1 x 10^^ eV and E'^Jy = 2.9 x 10 ^ GeV cm ^ s ^ sr ^, close 
to the upper limit "^-^-cascade" in Fig. [H This flux is detectable by JEM-EUSO, LOFAR 
and SKA. The cascade energy density does not exceed ou^"^. 

For cosmic F-and D-strings, the factor alone can increase the flux and Ucas to and 
above "E~^ cascade limit", keeping -Emin unchanged. 



H. EHE neutrino detection 



The calculated EHE neutrino fluxes are shown in Fig. [H together with the existing upper 



limits from ANITA 08 |43l, ANITA-II |43l, RICE[45[, and with the sensitivity of the proposed 



experiments - JEM-EUSO [37| (to be launched in 2015), LOFAR |39| and SKA |4q. The 
predicted fluxes are presented for gravity-dominated regime with z^, > z,^. As is discussed in 
Section IIV Gl in moduli dominated regime the neutrino fluxes can be also detectable. The 
characteristic feature of all models is very high minimal energy of neutrinos in spectrum 
i^min- Because of this, neutrino fluxes in some models are undetectable by JEM-EUSO. 

In all cases neutrino fluxes and cascade energy density are scaled by factor m^^^^. 
The quantity ^^"^ alone can be increased by factor 10-30. In our models neutrinos 

are produced at large redshifts, while ijJ^^ = 5.8 x 10~^ eV/cm^ is obtained mainly due to 
observation of 100 GeV photons [49*1, which have a local origin. Increasing the value of 
is limited since -Emin is usually proportional to and these high energy neutrinos become 
undetectable by JEM-EUSO. The quantity can be easily increased by factor 10^ for 
cosmic F- and D-strings. This factor is limited by Wcas only, -Emin does not depend on p~^. 

The characteristic feature of our model is the production of EHE neutrinos due to tremen- 
dously large Lorentz factors. The neutrino-induced gigantic showers in the air or lunar 
regolith produce a strong signal reliably detectable in optical and radio emissions. The sig- 
nature of EHE neutrinos with energies above 10^^ — 10^^ eV is given by this energy scale, 
which is unaccessible for astrophysical sources. 

This model has another signature, already discussed in |^]: the simultaneous appearance 
of a few showers in the field of view of a detector. It is due to neutrino propagation in 
the form of a very narrow jet and to the absence of time-delay in neutrino arrival, because 
of the tremendously large neutrino Lorentz factors = Ey/rriy. Compared to the case of 
superconducting strings in Q], the rate of multiple showers is strongly suppressed by higher 
energies of neutrinos, being partly compensated by the greater target mass of gigantic radio 
detectors. We present some brief estimates below. 

The number of detected neutrinos from a jet in the detector target is given by 

Nt'\> E) = ^-^(>-^) M,e,F.(> E), (77) 

where (Jun{> E) is the neutrino- nucleon cross-section for neutrinos with energies greater 
than E, Mdet is the target mass for the largest radio-detectors, and F„{> E) is the fluence 
of neutrinos with energy greater than E m. a. jet from a source at redshift z. The fluence is 
calculated as 

F.(> E) = j dfcAr^(fc) ^-^^'^^y-'°\ (78) 

where E^\-a is the minimum neutrino energy for a source at redshift z. In numerical estimates 
we shall use a^M ~ 1 x 10~^^ cm^ and Mdet ~ 10^^ g . 
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FIG. 1: Calculated neutrino fluxes compared with existing upper limits, ANITA-II, RICE, ANITA 
08, and with sensitivities of projects, JEM-EUSO, LOFAR and SKA. Line "E~^ cascade" presents 
the upper limit for cosmogenic neutrinos [49]. Curve "theor.3" gives the predicted neutrino flux 
for gravity-dominated regime with "normahzing" set of parameters p = 1 (ordinary strings), and 
Ol^ = ms = z^/z,y = 1. This flux is detectable only by SKA. Curve "theor.2" corresponds to 
superstring-motivated example in gravity-dominated regime (see Section HVFp withp = 1, z^/zi, = 
1, 07 = 3 and ms = 4. This flux is detectable by SKA and LOFAR. Increasing further the 
product ayrrig the flux can be made detectable by JEM-EUSO, though increasing of 07 should 
be limited, because of increasing of -Emin due to this factor. The curve "theor.l" with parameters 
p = z^/zjy = 1, ai = 2 and = 0.1 gives the flux detectable by all three future instruments 
SKA, LOFAR and JEM-EUSO. In case of moduli-dominated regime with parameters considered in 
Section HV Gl as aj = 2, = 0.1, and Z:^ = 100 neutrino flux is observable by all three detectors, 
SKA, LOFAR and JEM-EUSO. 



For the gravity-dominated regime, standard calculations give the following expression for 
the neutrino fluence from a loop of length Lmin at redshift z: 

^ 37r25/3 6, (tom)V2 [(1 + ^)i/2 _ i]2(i + ^)9/4 ^^.^'^p" 

where for z < 1 -Emin = 4.5 x 10^^m5^^yU_2o£ eV . Estimating fluence for the illustrative case 
parameters a = 3 x 10^, m = 4 x 10^ GeV, /i_2o = 10, and using 2; ~ 1 and -Emin ~ 10^'^ GeV, 
we obtain F,^ ~ 3 x 10^^^ cm~^. The number of neutrinos detected per burst is calculated as 
]\fdet 2 X 10~^n^ig. The fact that N^'^^ <^ 1 shows that practically all detected neutrinos 
are single. This is mostly due to the large value of -Emin in the gravity-dominated regime. 

The situation is different for the moduli-dominated regime ai < z < z^. The fluence 
from a loop with L ~ Lmin at redshift 2* is given by 
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For the model with a-j = 10, = 1 and z^, = 40, which yields the e-m cascade flux at the 
energy density bound oj'^f^, and assuming the neutrino energy E = 1 x 10^° eV, we obtain 
numerically Fu{> E) = 6.2 x 10~^^ cm~^, and the mean number of neutrinos detected in a 
burst is iV^^* ~ 0.4. The Poisson probability to detect simultaneously n = 2 neutrinos at 
average N = 0.4 is 0.054. 



V. CONCLUSIONS 

Production of high energy particles is a natural feature and one of the signatures of 
topological defects, including cosmic strings. It provides a method of searching for, e.g., 
cosmic strings, which complements other methods, based on the gravitational effects of 
strings, such as structure formation, CMB data, gravitational radiation, gravitational lensing 
and others. The strongest current bounds on strings with a symmetry breaking energy scale 
V = y/J^ is given by Gfi < 10~^ due to lensing effect [sil] and Gfi < 4 x 10~^, due to the 
millisecond pulsar observations [sH. With more advanced gravitational wave detectors, the 



bound is expected to improve to Gfi ~ 10-^2 [35|^ |53|. On the other hand, UHE particles 
can be detected from strings with G/i values as small as ~ 10~^°. 

Cosmic strings can arise from a symmetry breaking phase transition in the early universe. 
Fundamental strings of superstring theory can also play the role of cosmic strings in some 
models. 

A characteristic feature of the string dynamics is the periodic appearance of cusps, where 
very large Lorentz factors are reached for brief periods of time. Particle emission from cusp 
segment results in extremely high particle energies. Astrophysical sources, at the present 
level of knowledge, cannot accelerate particles to energies above lO^^ - lO^^ eV, with the 
maximum neutrino energy an order of magnitude lower. Thus, detection of neutrinos with 
E > 10^^ eV would be a signature of top-down models, with the string-cusp model as a 
plausible candidate. An additional signature of this model is neutrino emission in the form 
of a narrow jet with simultaneous detection of two or more neutrinos possible in the fleld 
of view of the detector. However, this possibility exists only for some values of model 
parameters and only for large detectors. 

EHE neutrino astronomy with E > 10^^ eV can probe high energy processes in the uni- 
verse up to red-shifts ~ 200. At the same time these neutrinos have large interaction 
cross-section with nucleons, a^N ~ 1 x 10~^^ cm^, favorable for detection. 

Moduli are produced near cusps of oscillating string loops, where the characteristic fre- 
quency of string motion exceeds the modulus mass m. These particles are emitted from 
a string segment with a Lorentz factor 7c ~ VmL/A, where L is the length of the loop. 
Moduli and the products of their decays move as a jet with an opening angle ~ 1c^- 
A modulus decays into two gluons, which initiate a quark-gluon cascade, which turns into 
hadrons, mostly pions, and then to neutrinos. We assume that the neutrino spectrum in 
the rest frame of the modulus is oc E~'^, where E^, is the neutrino energy in this frame. We 
adopted this spectrum in order to simplify the analysis and to obtain transparent analytic 
results. In fact the neutrino spectrum is not power-law, and has a flattening at a low energy 
e™^'^, which we call the minimum energy. This spectrum is Lorentz boosted, as described 
in Section IIII CI and the boosted spectrum also has a flattening at energy Emm ~ 7ce™°- 
This energy can be considered as an effective low-energy cutoff of the spectrum, because at 
E < -Emin the detectability of neutrinos is noticeably reduced. 

Given the length distribution of loops at epoch z and the spectrum of moduli emitted 
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from a loop of a given length, it is possible to calculate the neutrino flux E'^Jy{E), the 
cascade energy density Wcas, which provides an upper limit on the neutrino flux (see Section 
HVEl) . and the cutoff energy -Emin in the neutrino spectrum. The results are given in terms 
of three free parameters = a/10'', = GeV and /i_2o = G'/i/10~^°. Another 

important parameter is the redshift of the neutrino horizon, z,^ ~ 200. 

The results of these calculations critically depend on the redshift ^ AOOaj^^^ml^^ fi^J2o- 
This redshift separates two regimes in the string evolution: the gravity-dominated regime 
at z < z^, when gravitational energy losses are dominant, and the moduli-dominated regime 
at z > z*, when modulus radiation energy losses dominate. 

Apart from this selective role, fixing the value of gives a constraint in the parameter 

space (^7, 777,5, fJ^-2o)y in the form z^^a^, m^, fi-20) = where z'^ is the fixed z^: value. This 

reduces the number of free parameters to two (at the fixed z^), which we choose as aj and 

7775. As a result all calculated fluxes E'^J^{E) and energy density Ucas scale (at fixed 2;*) as 
-1 2 -1/2 

In our calculations the fixed value of z^ plays the role of the most important parameter, 
which determines what we call "the model". The high z^ models, defined as z^, > Zy corre- 
spond to the gravity-dominated regime in the whole allowed interval of redshifts (0, z^). The 
low z^ models, defined as < z^, are characterised by two regimes: the gravity dominated 
one at < z < and the moduli-dominated one at z^ < z < Zy, with approximately equal 
neutrino fluxes. The total neutrino flux is increasing with growth of z^. until it reaches Zj,, 
and models with z^ > Zy give the largest flux. This growth of flux with z^, is accompanied 
by a decrease of -Emin, see Eq. f l^ . which is favorable for neutrino detection by JEM-EUSO. 

Three theoretical curves in Fig. [T] illustrate different cases of neutrino detectability for 
the gravity-dominated regime at z^~^ Zy. The predicted flux is detectable by SKA only in 
case of "normalizing" parameter set p = a-j = = z^,/zy = 1 (curve "theor.3"). The flux 
is detectable by LOFAR and SKA in case of p = z^/zy = 1, = 3 and 7^75 = 4 (curve 
"theor.2"). The flux is delectable by all three detectors JEM-EUSO, LOFAR and SKA if 
p = z^jzy = 1, a-j = 2 and 7775 = 0.1 (curve "theor.l"). 

We considered above "ordinary" strings with p = \. For cosmic superstrings with re- 
connection probability p < 1, the neutrino flux increases by a factor p~^ without increasing 
Emini and is detectable for a wider range of model parameters. 

The cascade upper limit on neutrino flux in Fig. [T] [E^'^ curve) is given for cosmogenic 
neutrinos from [49!]. It is based on maximum energy density allowed by Fermi data 



48l |. For the the considered model with the dominant contribution from large redshifts this 
upper limit is higher because of 100 GeV-neutrino absorption. Therefore, allowed neutrino 
fluxes can be further increased. 

The remarkable feature of moduli-produced strings is strong beaming. The Lorentz-factor 
at z = Zy in gravity-dominated regime and at z^, in moduli-dominated regime is F ~ 10^^. 
The corresponding angle of a beam is ~ F"^ ~ 10~^^. Neutrinos with this tremendous 
energies arrive simultaneously at a detector and can produce simultaneously two or more 
showers in the field of view of detector. The estimates made in Section IIVHI shows that 
such events are rare due to very high energies of neutrinos. 
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APPENDIX 

In this section we shall derive the modulus radiation spectrum in more detail. In a 
flat background, i.e., g^u = Vfif = diag{— 1,1,1,1), the equation of motion for the string 
worldsheet X^{a"-) is 

da iV^Y'X^,) = 0. (81) 
Using the conformal gauge and = r, cx^ = a one obtains 

X/^ - X"^ = 0, (82) 

and the gauge conditions are 

X-X' = 0, (83) 

X^ + X'^ = 1. (84) 

In this gauge, the worldsheet coordinate r can be identified with the Minkowski time co- 
ordinate t. The solution for fl82l) can be written in terms of the right moving and the left 
moving waves as 

X(a,r) = l[X+(a+) + X_(a_)], (85) 

where the lightcone coordinates are defined as 0"+ = a + r, (T_ = a — r. The corresponding 
gauge conditions are X'+^ = X'_ = 1, where primes denote derivatives with respect to the 
lightcone coordinates. 

Using the lightcone coordinates, Eq.Q can be written in the form 

T(k,Un) = -- r da+ f da_ (l + X'^-X'_) e5[(-"-+-k-x+)-(-„<x_+k.x_)]_ ^gg^ 

L J-L J-L 

Since we shall be mainly interested in modulus bursts from cusps, we use the expansion 
of string worldsheet about a cusp, which we take to be at o"+ = (T_ = 0. The functions in 
the integrand of flS^ can be calculated from the expansions as 

l + XVX'_^-^a+a_, (87) 

and ^ 

k-X±^fc(^±a±T^4^, (88) 

where s is an 0(1) parameter which depends on the loop trajectory and k is assumed to be 
in the direction of the string velocity at the cusp. 
Eq. (!86l) can now be separated into two integrals as 

T(k,u;„) = ^^/+/-, (89) 
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where 

r-L 



/± = y da± a± e^n~"±"*"^"±J . (90) 
After a change of variables, we obtain the integral 

lAu) = ^ - l) £ci.xe-i4-^^1, (91) 

where 

and we have approximated the upper and lower limits of integration as ±00. The real part 
of the integral is zero since it is an odd function of x. The imaginary part can be expressed 
in terms of the modified Bessel function, 



Then, f l89|) can be calculated as 

T(k,a;„) = ^(^-l)V|/3(..), (94) 
and the power spectrum for the moduli radiation (jl]) from a cusp is 

The asymptotic form of the power spectrum for k » m and Un ~ k, i.e., u << 1, is^ 

"^^^ ta's'Gfi'n~'/\ (96) 



where F ~ 1. This is the same as the power spectrum for gravitons, except that for gravitons 
there is no additional coupling constant a and the numerical coefficient is somewhat different. 
The average rate of moduli radiation per solid angle is 

n " 

The sum over n can be converted into an integral over k by using the relation ujn = = 

y.=^l^^^- (98) 

n * 

Here we only consider the modulus bursts which have very large Lorentz factors, thus we 
keep the leading order term in the limit k » m. In this limit, ( 192|) becomes u ~ By 



3 When u « 1, K^{u) w (1) 
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substituting (195|) into ( 197|) . using ( p8|) and also by making a change of variable -u = -^^^ 
we obtain ^ ^ ^ 

rfiV ~ " K^/^{u)u^dudn. (99) 

The function K2/z{'u) dies out exponentially at m > 1. Hence, the main contribution to 
the rate comes from the region m < 1 which corresponds to 



k > kmin = kc^ -mVmL. (100) 

For A; >> k^nin, Eq- dHl]) gives 

dN ~ a^s^Gix^L^'^k'^/^dkdVt. (101) 

From f llOip . the number of moduli emitted from a cusp in a single burst, into solid angle 
dVL, having momentum between (/c, k + dk) can be estimated as 

dN ~ LdN ~ ah'^Gfi^L'^/^k-^/^dkdQ. (102) 

Here, we assumed one cusp event per oscillation period of a loop. 

Moduli are emitted into a narrow opening angle around the direction of the string velocity 
V at the cusp. The spectral expansion fllOll) has been calculated for moduli emitted in the 
direction of v. For moduli emitted at a small angle relative to v, Eq. fllOip still applies, 
but now the spectrum is cut off at kmax ~ l/Ld'^. In other words, the opening angle for the 
emission of particles with momenta ^k is 

~ {kL)-'/\ (103) 

Integration over Q in (llOip gives a factor ~ i?^, 

dN ~ ah'^Gfi^L-^/^k-'^/^dk. (104) 

The dominant contribution to the modulus emission comes from k ~ fcmm; and the total 
emission rate is 

The corresponding opening angle is 

^ ^ m/kmin- (106) 

The total power of the modulus radiation can be similarly calculated as 

P ~ a^.^G/.^L-V^A:-:f . (107) 
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